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ABSTRACT 

Stellar population synthesis (SPS) models are used to infer many galactic properties including star 
formation histories, metallicities, and stellar and dust masses. However, most SPS models neglect the 
effect of circumstellar dust shells around evolved stars and it is unclear to what extent they impact 
the analysis of SEDs. To overcome this shortcoming we have created a ne w set of circumstellar 
dust models, using the radiative transfer code DUSTY (Ivezic et al. 19991, for asymptotic giant 
branch (AGB) stars and incorporated them into the Flexible Stellar Population Synthesis code. The 
circumstellar dust models provide a good fit to individual AGB stars as well as the IR color-magnitude 
diagrams of the Large and Small Magellanic Clouds. IR luminosity functions from the Large and Small 
Magellanic Clouds are not well-fit by the 2008 Padova isochrones when coupled to our circumstellar 
dust models, and so we adjusted the lifetimes of AGB stars in the models to provide a match to the 
data. We show, in agreement with previous work, that circumstellar dust from AGB stars can make a 
significant contribution to the IR (> d/iw) emission from galaxies that contain relatively little diffuse 
dust, including low-metallicity and/or non-star forming galaxies. Our models provide a good fit to the 
mid-IR spectra of early-type galaxies. Circumstellar dust around AGB stars appears to have a small 
effect on the IR SEDs of metal-rich star-forming galaxies (i.e., when Ay ^ 0.1). Stellar population 
models that include circumstellar dust will be needed to accurately interpret data from the James 
Webb Space Telescope (JWST) and other IR facilities. 

Subject headings: stars: AGB and post-AGB - infrared: galaxies — galaxies: stellar content 


1. INTRODUCTION 

The physical structure, past history, and current prop¬ 
erties of a galaxy all go into shaping its observed spec¬ 
tral energy distribution (SED). As such, SEDs are pow¬ 
erful sources of information for unresolved galaxies and 
have long been used to discern the underlying physical 
properties of galaxies begi nnin g with the work of|Tinsley 
H972|),|Searle et al.|(|1973l), and Larson & Tinsley]p978 ). 
These studies pioneered the method of creating synthetic 
galactic spectra through the sum of the spectra of the 
stars hosted by the galaxy that has since become known 
as stellar population synthesis (SPS). By fitting galaxy 
SEDs with SPS models, properties of the galaxy such as 
the star formation rate, total mass in stars, metallicity, 
dust content, and the star formation history can be esti¬ 
mated. For details on the inner workings of SPS models 
and their b roader impact we refe r the reader to the re cent 
reviews by Walcher et al. (2011) and Gonroy|(2013). 

While there are exciting possibilities for the intorma- 
tion we are able to obtain through SPS model fitting, as 
Gharles Babbage noted when first introducing the me¬ 
chanical computer, we can only expect the right answers 
if we provide the ri ght input (or, c olloquially, “Garbage 
in, Garbage Out”; Babbage 1864). The limitations of 
inputs to SPS models are well-known, far ranging, and 
much discussed. Such limitations include incomplete 
isochrone tables and stellar libraries, poorly understood 
stellar evolution, and uncertainties in the initial mass 
function (IMF). In this paper, we focus on circumstellar 
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dust around asymptotic giant brach (AGB) stars, and 
the extent to which it affects the integrated light of stel¬ 
lar populations. 

The AGB phase is the last phase of stellar evolution 
in intermediate to low mass stars (~ 0.1 — 8Mq) in 
which significant nuclear burning takes place. During 
this phase stars eject their envelopes (with mass loss 
rates up to lO“'^M 0 yr“^) and evolve toward the white 
dwarf cooling sequence. T he material ar ound the star is 
observed to be dust rich (Bedijn 1987). This phase is 
notoriously difficult to model owing to the difficulty in 
modeling (in ID) stellar mass loss, convection, and mix¬ 
ing processes in the stellar interior. AGB stars are very 
luminous and can contribute many tens of percent to 
the integrated light of ste llar populations dMelbourne fc 


Boyer|20T3l Gonr^|2013[ [Melbourne et al.| 2012 [ [Kelson 


& Holden||2010|). 

I'he helium-shell flashes that occur during the AGB 
phase trigger the third dredge-up (TDU) process that 
brings material from the interior of the star to the sur¬ 
face. The surfaces of AGB stars begin as oxygen-rich 
(G/0 < 1) as a result of the primordial material from 
which the star forms, but the dredge-up process brings 
up carbon to the surface and in certain situations results 
in a carbon-rich envelope (C/0 > 1). The creation of 
Carbon stars depends sensitively on the various physical 
processes shaping the evolution of AGB stars. 

Despite these uncertainties, significant progress has 
been made in modeling AGB stars and their dust en¬ 
velopes and incorporating those models in SPS models. 
On the stellar evolution side, AGB models are becom¬ 
ing increasingly realisti c and constrained by observations 
(see Marigo et al.|2008 Girardi et al.|2010 Marigo et al. 
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Figure 1. Ratio of the flux including circumstellar dust to the initial input spectra for various values of the dust optical depth at l{im, 
tagb- The oxygen-rich models (left) are all for Teff = 2000 K while the carbon-rich models (right) are all for Tgff = 2400 K. The dominant 
feature in the oxygen-rich grid is the silicate feature at 10/im, seen in emission at moderate tagb values and in absorption at high tagb 
values. The silicon carbide feature at llfim is the dominant feature in the carbon-rich grid which also becomes more prominent for larger 
values of tagb • 


2013 

Cassara et al.||2013 

Rosenfield et al. 

2014 


have also been advances in understanding the gram prop¬ 
erties of dust arou nd AGB stars from detailed observa¬ 
tional studies (e.g. Suh 1999, 2000[ 2002 Groenewegen 


2006 


2011 


[Groenewegen et al.||20b9| |2012[ rinivasan et al. 

Sargent et al.||20lTy 


The connection between stellar evolution models, 
which make predictions for the photospheric properties 
of stars, and the associated dusty circumstellar envelopes 
is complex and highly uncertain. Nonetheless, several 
efforts have been made to integrate circumstel lar dust 
models into stellar isochrones and SPS models. IBressanI 
et al. (1998) adopted various empirical (or empirically- 
motivated; see e.g., jVassiliadis fc Wo^ |1993 Habing 
et al. 1994) relations between basic stellar parameters 


(i.e., mass, luminosity, and radius) and the mass-loss 
rate, pulsation period, velocity of the ejected material, 
and the optical d epth of the circumste llar envelope. The 
basic approach ofjBressan et al.|(|19 9 8|) has been updated 
and rehned by Piovan et al. l2008|), Gonzalez-Lopezlira 


et al. (2010) and Gassara et al. 


The influence ot A GB dust on the SEDs of gala xies is 


not well understood . Kelson & Holden (2010) and Chis- 
ari & Kelson (|2012 ) suggested that dust from AGB stars 


may be the dommant contributor to the mid-IR light 
in sta r-forming galaxies. However, Melbourne fc Boyerj 
(2013) argued that the former studies overstated the ef- 
te ct of AGB stars o n galaxies. In addition, the models 
of Silva et al. (1998) suggest that AGB dust plays a mi- 
nor role in the SHUs of actively star-forming (and star- 
burst) galaxies. In old stellar systems, it has been argued 
that the mi d-IR shows evidence for dust around evo lved 
stars (e.g., Athey et al. 2002 Martini et al. 2013). It 
is essential to understand these issues because the mid- 
IR is frequently used as a p roxy for the star formation 
rate (SFR) in galaxies (e.g., Kennicutt & Evans 2012). 
However, the connection between mid-IR tlux and bPR 


may not be quite as simple as commonly assumed, es¬ 
pecially for older stellar populations where heating of 
dust by older stella r populations (e.g., Salim et ar]|2009 


Utomo et al. 2014), and emis sion due to circums tellar 
dust around evolved stars (e.g. Martini et al. 2013) may 
play an important role. At z ~ 2, the Spitzer 24/im 
flux probes the rest frame ^8^m, and so understanding 
the contributions to the observed SEDs at ~10-30/im is 
critical for deriving reliable SFRs. In addition, it has 
been suggested that the mid-IR emission associated with 
AGB stars could be employed a s a stellar population age 
diagnostic ( ]Bressan et al.|[l9^ ). 

In this work we create rnodels for the dusty circumstel¬ 
lar envelopes around AGB stars and include those mod- 
els in the Flexible Stellar Population Syn thesis (FSPS 
Gonroy et al.||2009| |Gonroy & Gunn||20i0|) models. We 
employ t he publicly avai lable radiative transfer code 
DUSTY ( Ivezic et al.|1999 ) to model the dusty envelopes 
and then an empirically motivated prescription to assign 
dust shells to AGB stars in isochrones. 

The rest of paper is organized as follows: Section 2 
details the modeling of circumstellar dust shells and how 
they are connected to isochrones. In Section 3 we test 
and calibrate the new dust models using IR data of the 
Large and Small Magellanic Clouds. In Section 4 we ex¬ 
plore the parameter space of the models to see in which 
regimes the AGB dust has a signihcant influence on the 
integrated SEDs of composite stellar populations and 
compare the models to a variety of extragalactic data. 
Finally, in Section 5 we discuss the limits and implica¬ 
tions of our results. 

2. INCLUDING CIRCUMSTELLAR AGB DUST IN STELLAR 

POPULATION SYNTHESIS MODELS 

2.1. Modeling Dusty Envelopes 


We use the radiative transfer c ode DUSTY (Ivezic 
& Elitzur 1997 Ivezic et al. 1999) to model dusty cir- 
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Table 1 

Summary of the parameters chosen as input for the 
DUSTY models 


Parameter 

Garbon-Rich 

Oxygen-Rich 

Photosphere Properties 

Photosphere Model 

Aringer 

BaSeL 

Tfr 

2400 - 4000 K 

2000 - 4000 K 

log g 

0.0 

0.0 

Metallicity 

Zo 

Zo 

Envelope Properties 

Rjn Temperature (Tc) 

1100 K 

700 K 

Density Profile 

r-2 

r-2 

Dust Grain Properties 

K 

3200 fj,m?g ^ 

3000 iJ,m?g ^ 

Qext 

0.1 

0.1 

Grain Size (a) 

0.1 /im 

0.1 /im 

Density (pd) 

2.26 g cm~^ 

2.5 g cm~^ 


cumstellar shells around AGB stars. DUSTY solves the 
ID radiative transfer equation for a source embedded in 
a dusty region assuming spherical symmetry. The in¬ 
put to DUSTY is straightforward — it requires an in¬ 
put spectrum, the dust properties (chemical composition, 
grain size distribution, the dust temperature at the in¬ 
ner boundary), a radial density profile, and the optical 
depth of the envelope at a reference wavelength. We al¬ 
tered the wavelength grid to be finer than the default 
(1968 wavelength points instead of 105). In this section 
we describe our choice of inputs and summarize those 
choices in Table [T] 

We created two separate sets of models, one each for 
the carbon-rich and oxygen-rich AGB stars. In addition 
to composition, the models are a function of the optical 
depth, tagb (herein this refers to optical depth at 
and the effective temperature of the star. Teg. The model 
DUSTY spectra are implemented in the stellar popula¬ 
tion synthesis code differentially, in that only the ratio of 
the output to input spectra are stored and used within 
the code. 

In Figure we show the behavior of the oxygen-rich 
and carbon-rich grids as a function of optical depth at 
constant effective temperature for the carbon-rich and 
oxygen-rich grids. In this figure, we can see the char¬ 
acteristic features of the dust grains used in each set of 
models. For example, the 11/rm feature caused by SiG 
in the carbon-rich grid (right) and the lO^m silicate fea¬ 
ture in the oxygen-rich grid (left). The lOfim silicate 
feature varies from an emission to an absorption feature 
as a function of tagb- 

The differential implementation of the models allows 
us to keep the stellar parameters beyond effective tem¬ 
perature, such as stellar metallicity and surface gravity, 
constant. In Figure]^ we show the sensitivity of the dif¬ 
ferential spectra to surface gravity (middle panel) and 
metallicity (bottom panel). As can been seen, variations 
of these parameters does not have a significant impact on 
the differential spectrum. As such, surface gravity is kept 
constant at log(g)= 0.0 and the metallicity at Z = Zq for 
both the carbon-rich and oxygen-rich differential spectra. 

For the oxygen-rich input stellar spectra we use the 





Figure 2. Ratio of the flux including circumstellar dust to the 
initial input spectra for different grain distributions (top), surface 
gravity values (middle), and metallicity (bottom) to demonstrate 
how sensitive the AGB dust models are to these values. The red 
line in each panel represents our fiducial input spectrum used in the 
grid with a chosen effective temperature of 3000 K, all models are 
for tags =0.1. The differential spectra are generally insensitive 
to varying the parameters of the input spectrum. 


BaSeL (Lejeune et al. 1997) spectral library with effec¬ 
tive teniperatures ranging from 2000 K t o 4000 K. For 


the ca rbon-rich models we use the Aringer (Aringer et al. 


2009) spectral library with effective temperatures rang- 
ing from 2400 K to 4000 K. Stars hotter than ~ 4000K 
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Table 2 

Dust condensation temperature (Tc) 


Carbon-Rich Oxygen-Rich 

This work 

1100 K 

700 K 

Cassara et al. (2013) 

1500 K 

1000 K 

Groenewegen et al. (2009a) 

900-1200 K 

900-1200 K 

Groenewegen et al. (2009b)^ 

1000-1200 K 800-1000 K 

Marigo et al. (2008) 

800-1500 K 

800-1500 K 

Groenewegen (2006) 

1000-1200 K 1000-1500 K 

Piovan et al. (2003) 

1000 K 

1000 K 

Lorenz-Martins & Pompeia (2001)“* 

- 

417 - 1011 K 

Suh (2000) 

1000 K 

- 

Suh (1999) 

- 

700/1000 K ** 

Bressan et al. (1998) 

1500 K 

1000 K 

David Papoular (1990) 

- 

500 K - 800 K 

Rowen-Robinson Sz Harris (1982) 

- 

500/1000 K 


Values wer e fitted as part of a model 


Suh 


literature but chooses 1000 K for their own models 

are not expected to contain significant circumstellar dust 
shells. 

A key characteristic of the dust is quantified through 
the value of the extinction coefficient, k. 


Tlf^TTd Qext 

Pd 


( 1 ) 


where Qext is the extinction efficiency factor, a is the 
grain size, pd is the internal grain density, and nd is the 
spatial number density of a particular grain species in a 
dust mixture. The final k for each grain mixture is the 
sum of the k values for each component in t hat mixture . 
We fix Qext = 0.1 throughout, motivated by Suh (|1999|), 
and the grain densities are adopted from |Bressan et al.| 
([T9M|) (See Table 0. 

I'hroughout this work we adopt a grain size distri¬ 
bution that is a delta function at O.l^m. This is a 
common, though rather simpli stic assumption (e.g., |Suh 


1999 2002 


Piovan et al.||2003). We tested the impact of 


adopting diilerent grain size distributions in the DUSTY 
models and found that we can obtain similar emergent 
SEDs for different grain size distributions by varying 
TAGB (see top panel of Figure [^. 

In our model the parameters k, the dust-to-gas ratio, 
and i?in depend on whether the star is oxygen-rich or 
carbon-rich. To compute k we need to assume a certain 
dust composition. It is not currently possible to compute 
K from first principles owing to various uncertainties such 
as the detailed properties and evolution of circumstellar 
dust shells and dust grain formation and destruction pro- 


going 

(e-A; 

Jones et al. |2012l Ventura et al.||2014 

Nanni 

et al. 

2U13 

bchneider ef al.||2U14 

Dell'Agh et al. 

2UU 

)- 


As ICassara et al. 


discuss m detail, we expect 
the dust composition to change as a function of tagb ■ In 
principle, t his would require u s to compute n and tagb 


iteratively (Piovan et al.|2003). However, AGB stars can 
have diverse features from one another and it is unlikely 
that any one set of uniform models will fit all AGB stars 
equally. Since we have no a priori model for how the grain 
properties should vary with stellar properties, we choose 
to adopt a relatively simple, observationally motivated 
scheme. 


Sargent et al. (20101) and Srinivasan et al. (2010) fit 
models to mufti-wavelength broadband photometry of 
AGB stars with grain composition as a free parameter of 
the model. For oxygen-rich stars, | Sargent et al.| (|2010[ 
found that a grain composition of 100% oxygen-deficieni 
silicates was sufficient to fit the data. However, as an 
oxygen-rich shell becomes more opt ically thic k we ex¬ 
pect the silicates to become colder (Suh 1999 and Suh 


2002). Therefore, in our grid we adopt a gram type that 


varies with tagb- In our grid, for oxygen-rich stars we 
use warm silicates for the dust composition up until tagb 
= 3 when we start to include cold silicates in the mix¬ 
ture. The grain densities, sizes, and Qext parameters are 
all kept constant between the warm and cold silicates. 
For this reason, the extinction coefficient (Equation 1) 
remains the same as the grid transitions from cold to 
warm silicates. 


(1999 1 mentions that both temperatures have been used in the 


For carbon-rich stars, Srinivasan et al. (2010) found 
that a mixture of amorphous carbon (amG) and 10% 
silicon carbide (SiG) was sufficient. These dust composi¬ 
tions wer e used for their entir e gr id of AGB dust models , 
GRAMS (Sargent et al.|2011 and Srinivasan et al.|2011|. 
Furthermore, we found that chaiiging the amG and SiG 
ratio made only a modest difference in our results so we 
decided to keep the grain mixture for carbon-rich stars 
constant at 90% amG and 10% SiC for all values of tagb- 
Despite the simplicity of our dust compositions we 
find that they are sufficient for modeling individual AGB 
stars. We demonstrate this in Figure where we show 
how the DUSTY models compare to the SEDs of well- 
studied AGB stars HV 5715, SSTISAGEIC J052206.92, 
LPV 28579, and GW Leo. For the first three stars the 
data is as pr e sented in Sargent et al. ([20l0 ) and Sn^ 


vasan et al. (2010) with UBVl (Zaritsky et al. 1997 


Magellanic Glo uds Photometric Survey), JHK ( |gkrut- 


skie et al.|2006 2 Micr on All Sky Survey), Spitzer IRAG 
and MIPS bandpasses (Meixner et al. 1200^ SAG E), and 


Spitzer IRS spectroscopy (|Kemper et ahU^Ibl SAGE- 
Spec). GW Leo is an extensively observed object. Here 
we fit optical to IR photo metry (as present ed in Groe- 
newegen et al. 2012) - gri (Ahn et al. 2012 Sloan Dfgy 


ital Sky Sur ve^ V RI (|Le Bertre||1987|), JHKLM (|Le 
Bertre 1992), IRAS (|Beichman et al.|1988), AKARI (Ra 
et al. 1120081 AKARIIRG Survey of the Large Magellanic 


2008 %VKARI IRG Survey of the Large Magellanic 
Gloud), a.nd SP IRE and PACS bandpasses (|Groene\^ 


|gen et alT]|2011[ Mass-Loss of Evolved StarS). 

HV 5715, SSTISAGElG J052206.92, LPV 28579, and 
GW Leo are all oxygen-rich A GB stars and LPV 28579 
is carbon-rich. Sargent et al. (2011) presented detailed 


models of HV 5715 and SS'LISAgEIC J052206.92, Srini¬ 


vasan et al. (12011 ) modeled LPV 28579, and 


gen et ai. 12012) modeled GW Leo. The models used in 


Groenewe- 


h igure [Spwere made using the default parameters of our 
model grid and only allowing and tagb to vary. The 
shape of the observed SEDs and the observed features are 
fit well by the models. This comparison provides a test 
of our model grid both as a function of Tgg, tagb, and 
whether the star is oxygen or carbon-rich. We empha¬ 
size however that these fits are not unique and there are 
many degeneracies between e.g., the shell density profile, 
Tagb, grain size distribution, etc. 

The temperature of the dust at the inner radius, i?in, of 
the shell, Tc, is an important source of uncertainty in the 
modeling. Tablesummarizes the different values for Tc 













































































































































































Circumstellar Dust Around Evolved Stars 


5 



A ifim) 




A ifim) A ifim) 



A ifim) 


Figure 3. By-eye fits to four AGB stars using radiative transfer models. The blue dashed line is the input stellar photosphere and the 
red line shows the effect of including a circumstellar dust shell. The normalization in each panel is arbitrary. Far left: M star HV5715 
with Tefj = 3500K and a grain composition off 100% warm silicates. Middle left: M star SSTISAGEIC J052206.92 with Tefr = 3500K and 
grain composition of 100% warm silicates. Middle right: M star GW Leo with Teff = 2000K; this star is very dust enshrouded and so we 
adopt a grain composition of 100% cold silicates. Different lines are plotted showing different values of the dust temperature at the inner 
boundary, Tc- We show Tc = TOOK, lOOOK, 1300K (see text for details). There is a general agreement of the models for different Tc values 
especially at the mid-IR wavelengths. Far right: Carbon star LPV 28579 with T^ff = 3200K. 


adopted in the literature and demonstrates the range of 
values considered by various autho rs. Our choice for Tr . 
was based on fitting the data from Martini et al. (20131 
(as presented in Section 4). 

We compared the dust models of individual stars for 
the different values to the stars shown in Figure to 
check that on an individual basis the values gave reason¬ 
able agreement to AGB stars. In Figure we compare 
the different models for the range of values to CW 
Leo. Since CW Leo is the most dust enshrouded AGB 
star of our sample any effects due to changes to the in¬ 
put parameters will be amplified compared to the other 
stars in the sample. These comparisons lead us to adopt 
Tc = 700 K for the oxygen-rich grid. We emphasize that 
Tc could in principle vary with stellar type (i.e., with 
Teff), metallicity, etc., and our choice of a constant value 
for Tc may introduce systematic uncertainties in the final 
model results. 

The grid and the code used to generate the grid are 
publicly availably 

2.2. Connecting Dusty Envelopes to Stellar Isochrones 

With grids of circumstellar AGB dust emission as a 
function of optical depth, Teff, and C/0 taken from the 
dustless isochrones, our goal now is to connect the grids 
to stellar isochrones by computing the value of tagb at 
each isochrone point. 

The optical depth, tagbj is the key quantity that con¬ 
nects the stellar parameters provided by the isochrones 
to the models of the dust shells. Computing tagb from 
stellar parameters has been discusse d extensively in the 
l iterat ure b eginning with the wor k of |Vassiliadis fc Wood, 
( 19931) and iHa b ing et ah] (ll994|). Subseque nt wo rk by 
Bressan et al.| (|l 998|), IPiovan et ^ (|2003|, and |Cas-' 
sara et al. (|2(J13 ) rehned this technique to incorporate 
dusty envelopes into their SPS codes. In this work, we 
largely follow these previous efforts in coupling circum¬ 
stellar dust models to stellar isochrones. We will be brief 
in our description of computing tagb as the former three 
papers provide extensive details of the derivations of the 
following equations. 


^ https://github.com/AlexaVillaume/AGBGrid 
commit: 7308b5c424268cl6da3e4aeee9aef0b5f0bfcaf6 


We start with an initial equation for tagb that is 
derived assuming spherical symmetry and by integrat¬ 
ing over the thickness of the shell while assuming that 
Rout ^ 7?i„, where i?out and i?in are the outer and inner 
radii of the dust envelope: 


tagb = 


SMk. 1 
dtrCexp R'm 


( 2 ) 


In|Vassiliadis & Wood| (|1993|) they used observationally 
estirnated mass-loss rates of Galactic Mira variables and 
OH/IR stars to empirically determine equations for M, 
pulsation period (P), and Vgxp' 


'^exp 

km 


-13.5-b 0.056 


P 

days ’ 


( 3 ) 


with an additional condition that Wexp be in the range of 
3 — 15 km s“^, and: 


P R M 

log--=-2.07-(- 1.94 log--0.9 log—. (4) 

Q3,ys -^O -^0 

The equations for M are a function of initial mass 
of the star and whether the star is i n a su per-wind 
phase, defined by Vassiliadis & Wood (1993) as stars 
with P > 500 days. When not in the super-wind phase 
the mass-loss rate for a star with an initial mass < 2.5Mq 
is described by a simple relation with pulsation period. 


log M =-11.4-b0.0123P, (5) 

where P is in days, and M is in Mq yr~^, and for stars 
with initial mass > 2.5Mq, 


log M = -11.4-b 0.0125 


P- 100 


M 


- 2.5 


■ ( 6 ) 


In the super-wind phase the mass-loss rate is described 
as, 

M=— A. ( 7 ) 

C Vexp 


where c and Vexp are in units of km s 


-1 
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Figure 4. Top Panels: Cumulative luminosity functions for AGB stars for different ages and metallicities. Bottom Panels: Value of tagb 
as a function of Lbob metallicity. Only stars with Lboi ^ IO^Lq are shown. Stars that have the highest tagb values, i.e. the 

most dust-enshrouded, have the greatest bolometric luminosity, but are rare. 

The inner radius, i?in, is derived by equating the stellar 
luminosity with the luminosity at the inner radius, re¬ 
membering that the inner radius is defined as the radius 
where the temperature is equal to the dust condensation 
temperature, Tc, 


Rin = 


L 


47rT4 


( 8 ) 


We adopt our values from Table to obtain a final 
relation for carbon-rich stars, 


Rin = 1.92 X 10^^ 
and for oxygen-rich stars, 

= 4.74 X 10^2 


L 

Lq 


L 

Lq 


cm. 


( 9 ) 


( 10 ) 


The dust-to-gas ratio, S, is obtained by inver ting the 
observed cor relation between it and r;exp found in Habing 
et al. (1994), 


<5 = 


<5agb^^Lp / R 


225 


VlO-iLc 


-0.06 


( 11 ) 


where we adopt i5agb = 0.01 for oxygen-rich stars ( Suh 
19991) and <5agb = 0.0025 for carbon-rich stars (jhilanco 


et al. 1998). Note that in the following we assume that 
d does not depend on metallicity. See Section 2.3 for 
details. 

With the above equations we can now explore the be¬ 
havior of tagb along the isochrones as a function of 
metallicity and age. This is shown in Figure where 
we show the number of AGB stars as a function of lumi¬ 
nosity (top) and their corresponding tagb values (bot¬ 
tom). In this figure we see that the knee of the lumi¬ 
nosity functions coincide with the upturn of the tagb 
values. This indicates that brighter stars are both rarer 
and have larger values of tagb- Only for these rare stars 
is there an appreciable level of dust obscuration. 

The models a.re now included in the FSPS population 
synthesis code ( Conroy fc Gunn|2010 ) as of v2.5 ^ FSPS 
computes tagb as per the equations above for identified 
AGB stars in the isochrones. With tagb, LeS, and the 
composition of the stellar envelope, FSPS interpolates 
within the grid of DUSTY models and grafts the inter¬ 
polated model onto the stellar SED. 

We briefly review here the salient characteristics of the 
FSPS stellar population synthesis model. FSPS takes as 
input a set of stella r isochrones (h e re we use the Padova 
models detailed in |Marigo et al.| (|2008|)), stellar spec¬ 
tral libraries (in our case the BaSeL library), allows the 
user to specify a stellar initial mass function (IMF), and 
outputs simple stellar populations (SSPs), and, if a star 

^ code.google.com/p/fsps 
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formation history is specified, composite stellar popu¬ 
lations. FSPS includes models for diffuse dust absorp¬ 
tion and emission. The model has been exte nsively cal¬ 
ibrated and tested against observations, see [Conroy k\ 
Gunn (2010) for details. 

When AGB dust is turned on in FSPS, tags is com¬ 
puted at each isochrone point using the method described 
above. Fundamentally, tags depends on the mass and 
age of the star. 

2.3. Metallicity Dependence of the Dust-to-Gas Ratio? 

To first order, the seed elements relevant for grain for¬ 
mation in oxygen-rich environments, such as O, Si, and 
Fe, will vary in proportion to the metallicity of the star. 
It therefore might seem reasonable for the dust-to-gas 
ratio to increase linearly with the metallicity of the star 
(e.g., Habing et al.||1994 |. 

However, recent models attempting to follow the for¬ 
mation and dust grains self-consistently as a function of 
AGB mass and metallicity have revealed a complicated 
relat ion between the dust-t o -gas ratio and metall icity 
(e.g. [Ferrarotti & Gail||2006|). [Schneider et al.|(|2014|) de- 
tailed the complicated factors that contribute to the dust 
production of AGB stars by comparing the observed dust 
production rates of AGB stars in the L arge and Small 
Magel lanic Clouds to theoretical models. Schneider et al.[ 
(2014) noted that efficiency of the Hot Bottom Burning 
(HBB), the TDU, and the mass of the star all contribute 
to the ef ficiency of the dust pr oduction for AGB stars. 
And, as Ventura et al. (2014) noted, the efficiency of 
these mechanisms are influenced by the modeling of con¬ 
vection in the star. Discrepancies between the theoretical 
dust yields a nd observed dust yields found in [Schneider 


et al. (2014) highlight that we do not have a solid pre¬ 


dictive model for dust formation in AGB winds. More¬ 
over, current models suggest that dust production could 
be a complex, non-monotonic function of metallicity. For 
these reasons we decided to adopt a simple approach and 
keep the dust-to-gas ratio fixed for all oxygen-rich stars, 
independently of metallicity. 

For carbon-rich stars the seed elements for the grains 
(carbon) are made in-situ in the star, so one would expect 
that the dust-to-gas ratio to have little dependence on the 
original metallicity of the star. Hence also for carbon-rich 
stars we keep the dust-to-gas ratio fixed, independent of 
metallicity. 

3. CALIBRATING STELLAR POPULATION MODELS THAT 
INCLUDE AGB DUST 

As an initial test of the new models, we compare 
them with the photometry of the Large Magellanic Cloud 
(LMC) and the Small Magellanic Cloud (SMC). The data 
is from the Sur vey of the Agents of Galaxy Evolution 
(SAGE) survey ([M eixner et al.[|2006 for LMC data and 
Gordon et aL|20lT for SIVIG). The SAGE survey is based 
on Spitzer IKAG and MIPS data, including photometry 
in [3.6], [4.5], [5.8], and [8.0] b ands. _ 

For the LMC data we follow jCassara et al. (2013) and 
select a region within tt square degrees of the center of 
the LMC (RA = 5h23m.5, DEC = -69,45’). For the 
[3.6] - [8.0] vs [8.0] color-magnitude diagram (CMD) for 
both galaxies we make a further cut on the dim, red 
background objects for clarity. In Figure [^ we include 
isochrones for both the LMC and SMC, with respec¬ 



Figure 5. Color magnitude diagrams of stars in the Small Magel¬ 
lanic Cloud (top) and the Large Magellanic Cloud (bottom) over¬ 
laid with SSP models of various ages. Foreground objects have been 
removed with a color-magnitude cut. The solid lines are single-age 
models that include AGB circumstellar dust while the dashed lines 
show the same models without AGB dust. Note that the luminous 
red spur in the data is only captured by the models that include 
AGB dust. 




Figure 6. Gomparison of the observed and model luminosity func¬ 
tions for the SMG and LMG. Model LFs were deri ved using the 
default Padova 2008 i sochrones JMar i go et al.|2008ll. the updated 
Padova isochro nes by iGirardi et al. ( Case A), isochrones 

as calibrated by|Conroy 66 tiunn|l|2Ul0ll, and isochrones calibrated 
for this work. All model predictions include circumstellar dust for 
AGB stars and have been convolved with estimated SFHs of the 
LMC and SMC. 

tive metallicities oi Z = 0.008 and Z = 0.004, from the 
Padova stellar models. 

In Figure there are six models in each panel: three 
ages, and for each age there are models both with and 
without AGB dust (solid and dashed lines, respectively). 
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Figure 7. SSP colors as a functi on of aee, comparing m odels 
from this work to those described in [Conroy fc Gunnl ||2010 i . This 


figure al so includes LMC st ar cluster data tromjlNoel et_aL 


120131 


(circles), jPessev et al.|(|2008|l (triangles), andlGonzalez et al.|l|2UU4 1 
(squares). By construction, the colors predicted by the c alibrations 
described in this work are very similar to the colors from jConroy fcj 
jGunnj poTojl. Notice that the very similar integrated light colors 
between the two models nonetheless produce very different LFs in 
Figure |6| 

Overall, the models with circumstellar dust cover the ob¬ 
served behavior of both the LMC and SMC, providing 
support to our method for assigni ng dusty envelopes to 
AGB stars in isochrones (see also Cassara et al. (20131 
who find similar agreement with a different set of ACB 
dust models). 

Comparing isochrones to data in CMD space hides po¬ 
tentially significant discrepancies in the relative numbers 
of stars in different evolutionary phases. We provide a 
closer comparison of our models to LMC and SMC lumi¬ 
nosity functions (LFs) in Figure]^ The observed cumu¬ 
lative LFs (thick bla ck line in Figure!^ are constructed 


from the catalogs of Boyer et al. (^2011), including only 


those stars classihed as any variety of O, C, or x-AGB. 
We have inspected images of the brightest objects to in¬ 
sure that they are not strongly affected by confusion or 
high backgrounds. 

We create model LFs by convolving model isochrones 


with measured SFHs and metallicity histories from Har- 


ris & Zaritskyl |2004|) and applying a cut of Teff < 4000 K 
to isolate the AgB stars. We checked that the convolu¬ 
tion of the SFHs with the stellar mass and integrated lu¬ 
minosities reproduces the observed integrated NIR light 
of the clouds and the to tal s tellar masses repor t ed by 


Harris & Zaritsky (2004) and Harris & Zaritsky (2009) 
to within a factor of two. We show four modefs in h'ig- 
urel^ three of which differ o nly in the input iso chrones: 
unaltered Padova isochrones (Marigo et al. 2008) (MG08, 
dashed blue line) , re-ca librated Padova isochrones from 
Gonroy & Gunn ( 2010[ ) (GGIO, dot-dashed green line). 


and a new normalization created specifically to provide a 
better fit to the observed LF (sold red line, described in 
detail below). The fourth is compu ted from the updated 


Padova isochrones as described in Girardi et al. 

(Gase A). We emphasize t hat this is an entirely different 


([Mo| 


model where we ad opt the Girardi et al . (2010) bolomet- 
ric corrections and|Marigo et al.|(|2008|) dust models. In 
all cases circumstellar dust around AGB stars has been 
included. 

Uncertainties in the SFHs are difficult to incorporate 
correctly due to the substantial covariance between the 
SFRs in adjacent temporal bins, but the uncertainties on 
the SFR in any given bin are typically less than a factor 


of two, and this gives an upper limit on the uncertainties 
in the predicted CLFs at any magnitude. The effects 
of photometric uncertainties on the predicted GLFs are 
negligible, though uncertainties substantially larger than 
reported may affect the bright end. 

T he predicted lumino sity functions deter mined from 


(2008) isochrones and the Gonroy & 
calibrations do not match the data well. 


the M arigo e t al. 

Gunnl (|2010 j 

Models oased on the Marigo et al. (2008) isochrones over- 
pre dict the number of AG B stars while models based on 
the Gonroy & Gunn (2010) normalization do not span the 
full extent of the rnagmtude range. Due to these short¬ 
coming we decided to introduce updated calibrations of 
the Padova isochr ones for the AGB phase, different from 
t hose presented in Gonroy & Gunn (2010). 

Conroy & Gunn (|2010l) chose to decrease the overall 
contribution of AGB stars to the integrated light by de¬ 
creasing log L and slightly modifying This modih- 
cation resulted in a better fit to the NIR colors of inter¬ 
mediate age star clusters in the LMC. For the new cal¬ 
ibration, we chose to keep log L unchanged and instead 
adjust the IMF weight to effectively reduce the lifetimes 
of AGB stars: 


weight = max(-k0.1, io-i o+(*-8 0)/2.5))^ (-^2) 

IMF^gight = weight x IMFweight, (13) 

where t is log age (yr). This change keeps t he inte¬ 
grated light pr edictions similar to the ordinal |Conroy| 
& Gunn|(2010|) normalization (see Figure [t] and below). 
i'his new calibration is the default re-normalization for 
Padova isochrones in FSPS as of v2.5. The calibration 
introduced in this work is an improvement from previous 
versions available in FSPS. 


The dusty isochrones from Girardi et al. (2010|) are 
comparable to the models presented in this work, 'ihere 
is a difference in the extent in magnitude that two mod¬ 
els span seen at 4.5/rm. At 8.0/iTO, however, the overall 
shape of the luminosity functions are similar. We com¬ 
pared all the dust schemes made available by the Padova 
group and found none that made an imp rovement the 
magnitude range in the Girardi et al. (2010) models. The 
model LMC luminosity function presented in this work is 
well-matched to the data with the exception of the bump 
seen at 7 < [8.0] < 9. 

However, while improved from previous versions, the 
model SMC luminosity function still significantly over¬ 
predicts the number of AGB stars. The source of this 
discrepancy is currently unclear but it could be a prob¬ 
lem with the treatment of the AGB evolution in the lower 
metallicity isochrones. Agreement between the models 
and data for the SMC in Figure [^indicates that we are 
likely not over-predicting the amount of AGB dust per 
Lboi but, again, the CMDs hides information regarding 
AGB lifetimes which are critical for the luminosity func¬ 
tions in Figu re Given that a similar over-prediction 
is seen in the Girardi et al. (2010) models indicates that 
this is an issue more tundamentai than the dust models. 

As a further check, we compare predicted AGB number 
counts for several galaxies from the ACS Nearby G alaxy 
Survey Treasury (ANGST [Dal canton et al.||2009| ) using 
the HST WFC3 F160W filter. iMelbourne et al.f(pOT2) 
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Figure 8. Model spectral energy distributions as a function of age comparing models with (solid lines) and without (dashed lines) 
circumstellar dust around AGB stars. Notice that circumstellar dust has a very large effect on the SED for A > 4/im. For the Padova 
stellar models that we employ, the effect of circumstellar dust around AGB stars is relatively modest at 10 Gyr. 


Table 3 

Comparison of AGB counts for ANGST galaxies 


Galaxy 

Observed 

MG08 

GIO 

This Work 

This Work 

Data 

DD078 

273 

858 

721 

724 

2.65 

DD082 

1046 

2556 

1122 

1709 

1.63 

IG2574-SGS 

1504 

2587 

1554 

1867 

1.24 

NGC4163 

640 

1425 

677 

945 

1.47 

UGC4305-1 

740 

1438 

1113 

1136 

1.54 

UGC4305-2 

721 

1532 

1063 

1428 

1.98 


compared the observed number of AGB stars for ANGST 


galaxies to the predictions by the ( 

Marigo et al. 

2008 

MG08) isochrones and the ( 

Girardi et al. 

20TU 

TITUj 


models that include the updated normalization of the 
AGB lifetimes in Table As expected, our model pre¬ 
dictions are lower than those of MG08. However, there 
are still cases where the discrepancy with observ ations 
is significant, as noted in Melbourne et al. (20121 when 


adopting the more recent Gil) isochrones. Uur mean 
model-data num ber ratio for AGB stars is 1.75, similar 
to the value from Melbourne et al. (20121 when using the 
GIO isochrones. 

Our new normalization of the Padova isochrones was 
chosen to provide a decent fit to both the LMG lumi¬ 
nosity function and the integrated colors of LMG star 
clusters. The latter is shown in Figure In this figure 
we compare both our new normalization and the normal¬ 


ization from G GIO to observations (See Figure 2, Conroy 
fc Gunn|2010 ) . The latter are obtained by stacking star 


clusters in bins of age to reduce the effects of sparse sam¬ 
pling of the AGB in individual clusters. By construction, 
the integrated colors are well-matched to the data. 

While we have calibrated our mode ls based on Fig- 
urepzl we note that Girardi et al. (2013) cautioned about 
biases introduced due to a “boosting” of the number 
of AGB stars right at the age of many LMG clusters. 
However, the good agreement between the observed and 


model LMG LF indicates that the effect of AGB boost¬ 
ing over a narrow range of ages is unlikely to have a 
significant effect on our results. 

We wish to emphasize that our renormalization of the 
Padova isochrones is completely ad hoc and is simply 
chosen to provide better fits to observations. Definitive 
conclusions regarding AGB lifetimes, etc., will require ac¬ 
tual modifications to the underl ying physics in the st ellar 
models, as for example done in Girardi et al. (2010). 


4. THE EFFECT OF AGB DUST ON INTEGRATED LIGHT 
4.1. General Trends 

In this section we first describe the general conditions 
in which circumstellar dust may affect the integrated 
light of simple and complex stellar populations (§4.1), 
then move to discussing two examples in which there 
is evidence for the impact of circumstellar dust (early- 
type galaxies in §4.2 and low-metallicity dwarf galaxies 
in §4.3), and one example in which circumstellar dust 
appears to play no important role (actively star-forming 
galaxies, in §4.4). 

With the new circumstellar dust models we can calcu¬ 
late revised SEDs and compare them to models without 
circumstellar dust. lu Figure 1^ we show model SEDs 
with (solid hues) and without (H^ashed lines) dust shells 
around AGB stars for ages of 0.2, 2, and 10 Gyr for 
Z = 0.2Zq and Z = The circumstellar dust sub¬ 
stantially increases the model flux red-ward of ~ 4 /xto. 
At old ages the influeuce of the AGB dust is sharply di- 
miuished, having little discernible impact on the SED. 
We note here that the models in this figure contain no 
diffuse dust associated with the interstellar medium (see 
below). 

The trends with age and metallicity in Figure can at 
least partially be understood by considering separately 
the effects of carbon-rich and oxygen-rich stars on the in¬ 
tegrated SED. We explore this in Figure which shows 
the fractional contribution of oxygen-rich and carbon- 
rich stars to Lboi (top panels) and the relative contribu- 
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log Age (yr) 


log Age (yr) 


log Age (yr) 


Figure 9. Relative contributions of carbon and oxygen-rich stars to the integrated light of single-age stellar populations. Top Panel: 
Fractional contribution of carbon and oxygen-rich stars to the total bolometric AGB luminosity. Bottom Panel: Relative contribution 
to total 24/i.m flux, comparing models with and without AGB dust. This flgure highlights both the age and metallicity dependence of 
carbon-rich star formation and also the disproportional influence of carbon-rich stars to the mid-IR flux when AGB dust is included in the 
models. 


tion of oxygen-rich and carbon-rich stars to the 24/xm 
flux (bottom panels). Note that we are using Padova 
stellar models here. 

We see the well-known fact that the appearance of 
carbon-rich stars is both metallicity and age dependent. 
In the top row of Figure we see that for the lowest 
metallicity model, Z = O.IZq, carbon-rich stars dom¬ 
inate Lboi for a longer period of time than the higher 
metallicity models. Moreover, carbon-rich stars do not 
exist at old ages, 9.3 < logt/yr < 10 for high metallici- 
ties. The age dependence is complicated, depending on 
processes such as the TDU and HBB. Detailing when and 
under what conditions carbon-rich stars are produced is 
outside the scope of this work but it is a much discus sed 


subject in the literatur e (e.g., Iben & Renzini 1983 
Criscienzo et al.||2013l and the references therein). 


Di 
Tt 

is important to understand because at certain ages and 
metallicities carbon-rich stars are dustier and therefore 
redder than oxygen-rich stars, as demonstrated in the 
lower panels of Figurej^ We see in Figurej^that carbon- 
rich stars have a disproportionate influence on the 24/rTO 
flux relative to their overall Lboi contribution. The lack 
of carbon-rich stars at old ages, combined with the gen¬ 
erally lower AGB luminosities at late times (see Figure]^ 
explains why circumstellar dust emission from AGB stars 
has a relatively modest effect on the integrated light SED 
at old ages. 


In Figure [T0| we show how AGB dust affects the NIR 
and mid-IR integrated light of composite stellar popula¬ 
tion models as a function of age for different metallicities 
(note that this refers to populations as a whole, rather 
than the metallicity of individual stars discussed in §2.1, 
which will have an effect on the SPS models), star forma¬ 
tion histories (SFHs; represented by exponential models 
with an e-folding time given by tsf), and diffuse dust 
content (the contribution of dust around young stars and 
dust in the ISM). In this figure four wavelength ranges 
are shown and in each panel we compute the ratio of 
luminosities between models with and without circum¬ 
stellar AGB dust. 

The diffuse dust model is a combination of the lGharlotl 
& Fall (2000| model for dust attenuation, with a power- 


law attenuation curve (with Tdig quoted at SSOOAand 
the power-law index, a = —0.7 unless stated other¬ 
wise) and the 
dust emission. 


(2007) model for diffuse 


Draine & Li 

i'he overall amount of diffuse dust emis¬ 
sion is determined by energy conservation. In this pa¬ 
per, models that include diffuse dust also include dust 
around young stars, the so-called “birth cloud compo- 
nent”, where The = 3rdiff. Note that the |Draine fc Li| 


(2007| dust models include both thermal emission from 
cold dust and PAH emission. The models are described 
by three parameters gpAH, Gmin, 7- Unless otherwise 
stated we adopt 3.5%, 1.0, 0.01 respectively for these 
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Figure 10. Effect of circumstellar dust on model fluxes as a function of metallicity (left column), star-formation history (middle column), 
diffuse dust content (right column), and wavelength (top to bottom). The model plotted in red is the default model, same across every 
panel, with Z = Zq, rsp = 1 Gyr, and r^iff = 0.0. The diffuse dust content significantly reduces the influence of AGB circumstellar dust 
at > Sfim but at 4.5^m the models are largely insensitive to the presence of diffuse dust. 
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Figure 11. A more detailed look at the effect of diffuse dust on 
the contribution of AGB dust at 4.5/im (top) and 24/im (bottom) 
for several ages. At 24/i.m there is dramatic decline in the influence 
of AGB dust for even moderate values of diffuse dust. However, 
flux at A.bfim is largely insensitive to the contributions of diffuse 
dust as the diffuse dust in our model is not hot and therefore does 
not radiate substantially at such short wavelengths. 


values. 

We see in all models that only after ^0.1 Gyr does 
AGB dust become influential, peaking at about 1 Gyr, 
with a significant reduction by the time the galaxy 
reaches 10 Gyr in age. The influence of AGB dust on 
the SED varies in interesting ways as a function of wave¬ 
length, metallicity, and SFH, but by far the most influ¬ 
ential parameter is the amount of diffuse dust. 

For galaxies where there is no diffuse dust the impact 
of the dust around AGB stars is significant and extends 
to wavelengths as short as ~ 4/im. However, even a 
modest amount of diffuse dust significantly diminishes 
the relative influence of AGB dust at > 8fim. At shorter 
wavelengths the effect of AGB dust dominates over that 
of diffuse dust, since in our model the diffuse dust com¬ 
ponent is generally not hot enough to emit significantly 
at < 8fim. We investigate this in more detail in Fig¬ 
ure 11 where we show the effect of AGB dust at d.S/rm 


and 24/rTO as a function of the amount of diffuse dust 
for several ages. At 4.5/rm the AGB flux is relatively 
insensitive to the diffuse dust content but at 24^m the 
diffuse dust greatly reduces the influence of the emission 
from the AGB dust shells. From ^ 1.0 — 10.0 Gyr the 
AGB dust has a significant influence of the mid-lR light 
only for low diffuse dust values and, consistent with the 
previous figures, drops off as the diffuse dust becomes 
more prominent. We therefore conclude that circumstel- 
lar dust around AGB stars will be important in regimes 
where the diffuse dust optical depth is roughly <0.1. 

We also investigate the effect of AGB dust on optical- 
NIR and NlR-mid-lR colors in Figure The difference 
is not large for V — Kg, but this color is commonly em¬ 
ployed in high-precision SED modeling of galaxies and so 
even a small change could potentially induce large differ¬ 
ences in derived parameters. The models that include 
AGB dust tend to be bluer in. V — Kg than models that 
do not for a range of metallicities and ages. This is be¬ 
cause AGB stars, with or without dust shells, contribute 
essentially zero flux in the V band and signihcant flux 
in Kg. Dusty AGB stars can be optically thick even in 
the Kg band, resulting in a lower flux and hence bluer 
V — Kg color. In the middle and bottom panels of Figure 
[T^ we see that the models with AGB dust are signif¬ 
icantly redder than models without for optical-mid-IR 
and mid-IR-mid-IR colors. 


4.2. Comparison To Early-Type Galaxies 

Early-type galaxies are potential candidate systems 
where the effect of AGB dust may be important, owing to 
their generally low diffuse dust content. They are known 
to host intermediate and old age populations, 1 < t < 10 


Tracer et aL||19981 |Thomas et al.||2005 IJimenez et al. 

2007| 

Trager et al.||2005| |Thomas et al.j|2005| 

Conroy 

et al. 

evider 

2005) 

exists 

gestin 

20141). While many early-type galaxies c 
ice for diffuse dust at low levels (e.g |Lau€ 
the origins of which are actively debatec 
galaxies with only upper limits in the FI 
g that they have very little diffuse dust (e 

0 show 

r et al. 
1, there 
R, sug- 
g. Jura 

et al. 111987 IKnapp et al.||1989| |Temi et al.||2009[ |Martini 


et al.||2013|. 


Athey et al. (2002) presented a mid-IR spectrum from 
the infrared Space Observatory (ISO) for the early-type 
galaxy NGG 1404. These authors found a clear excess in 
the data at 8—14 /xto compared to a blackbody extrapola¬ 
tion. They attributed this excess to 0-rich AGB features. 
In Figure 13 we show their data for NGG 1404 compared 
to our models with no dust (blue line) and three models 
with AGB dust for different ages (red lines). All models 
assume the galaxy has a metallicity of ^0. 

First, we see that by including AGB dust in the model 
the fits to the data are greatly improved at 8 — 12p,m. 
Furthermore, we see from the three models with AGB 
dust the mid-IR ca n be a powerful age di s criminator, as 
alread y disc ussed in Bregman et al. (2006), Bressan et al. 


(19981, and Athey et al. (|2002|), for simple systems with 
only one burst of star-formation. We find that a model 
with AGB dust and an age of ~7.5 Gyr provides the 
best fit to the 10pm feature. This is in good agre ement 
with the 9 Gyr estimated by Bregman et al. { 20061), who 


used the AGB dust models from Fiovan et ah 12003) to 
estimate galaxy ages based on niid-lR SFDs. This age 
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Figure 12. Effect of AGB dust on optical and NIR colors. The 
effect on mid-IR colors is large and very time- and metallicity- 
dependent. The contribution of circumstellar AGB dust to V—Ks 
is small compared to the other colors, but V—Kg is frequently used 
in high-precision SED modeling of galaxies and, with the typically 
small measurement uncertainties of optical-NIR colors, the rela¬ 
tively small effect of circumstellar dust on V—Kg may in fact have 
a significant effect on parameters derived from SED modeling. 


is also in agreement with an optical spectroscopically- 
estimated age for this galaxy of 7 Gyr (Conroy & van 
Dokkum, in prep). 

Interpretation of the mid-IR data for NGC 1404 is com¬ 
plicated by the fact that this galaxy is apparently not 



Figure 13. ISO spectroscopy for the early-type galaxy NGG 1404 
(points with errors), compared to models with (red) and without 
(blue) AGB dust. The model with AGB dust gives a much better 
fit to the data at > Sfim. Moreover, it is clear that the mid-IR is 
very sensitive to age even over the interval 6 — 10 Gyr. 


entirely devoid of diffu se dust. It is detec ted by Herschel 
at 70^m and 160/rm (Temi et al. 20091, indicating the 
presence of cold dust. Howev er, we attempted to fit the 
data from Athey et al. (2002) with only diffuse dust (no 
AGB dust) and found that such a model could not simul¬ 
taneously fit the mid-IR and FIR data (see also Figure 


15). Therefore, we conclude that circumstellar AGB dust 


appears to be necessary to fit the mid-IR data for this 
g alaxy. 


Martini et al. (2013) recently published mid-IR pho- 
tometry of ea rly-type galaxie s froni Spitzer's IRAC and 
MIPS bands. Martini et al. (2013) divided the sample 
based on whether there were diffuse dust lanes de tected 
in HST imag i ng an d detected at 70 and 160 fim. |Ram-| 


pazzo et al. (2013) also published an atlas of mid-IR 
spectroscopy from Spitzer-IRS of early-type galaxies that 
they sort into classifications based on the contribution of 
old stellar populations to the mid-IR continuum, where 
class 0 are galaxies where the old stellar populations are 
most signifi can t. 

In Figurej^we compare the data for galaxies that did 
n ot have detected dust lan es and were classified as class 0 
in|Rampazzo et al. (2013) to models that include (solid- 
red) and do not include (dashed-blue) AGB dust. All 
models plotted were made for a solar metallicity, single¬ 
age stellar population. The data were normalized to the 
flux at 5.72/rm, which should eliminate aperture size is¬ 
sues between IRAC and IRS. The ages of the galaxies 
were estimated by hand to provide a satisfactory fit to the 
3 — 20p,m data. The model that includes AGB dust fits 
the 3.6 — S.Ofim data just as well as the model that does 
not, only for > 8.0/im are there significant differences 
between models with and without circumstellar dust (for 
these ages). The models provide a good fit to the mid- 
IR spectra from Spitzer-IRS. In particular, the overall 
shape, from 5 — 20/im, is reproduced well and the rel¬ 
ative strength and positions of the broad 10 and 18pm 
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Figure 14. Spitzer IRAC and MIPS photometry and IRS spectro scopy of a sample of e arly-type galaxies co mpared to Z=Z0 models with 
(red) and without (blue-dashed) AGB dust. Photometry are from|Martini et aL|(|2013|l and spectra are from|Rampazzo et al.|||2013|l. All 
data were normalized to 5.72 pm, the IRAC 3 central wayelength. 1 he errors in the photometry are smaller than the points, i tie ages 
were determined through by-eye fits. The models that include AGB dust generally fit the data well, with slight oyer predictions of the the 
10pm emission feature seen in all the spectra. 


dust features are generally well matched. 

In Figure [T5| we compare the observations of NGC 2300 
with both the AGB dust model and two models that 
includes only diffuse dust. The fundamental difference 
between the latter two models is the adopted value of 
9PAH- In one case we adopt the fiducial value while in 
the other we adopt an unreasonably high value of 10% 
0 We added as much diffuse dust into the models as 
possible without exceeding the constraints of the FIR 
flux at 70/iTO and I60/im. It is clear that the model 
with only diffuse dust cannot reproduce the 8 — 24/rTO 
features seen in early-type galaxies with influential old 
stellar populations. Varying the diffuse dust parameters 
qpAH and Umin does not qualitatively improve the fits. 

The slight over prediction of the strength of the 10fj,m 

® This is a n extrapolation of th elDraine &: Lilll2007ll models, the 
highest value [Praine Li|||2007| compute is 4.&b%. 


emission feature is likely due to our treatment in the dust 
grain properties. We reiterate that our grain scheme is 
relatively simplistic and there is an abundance of evi¬ 
dence that the actual grain mixtures in AGB stars, and 
their dependence on the properties of the star, are more 
complica ted than represented in the models pr esented 
here (e.g. Jones et al.|20I2 McDonald et al. 12010 1. Other 
groups producing siniilar dust shell models have incorpo- 
rated mo re detailed approache s to the grain schemes. For 
example, Cassara et al. (2013) detail a sophisticated ap¬ 
proach to the grain mixtures that account s for changes 
in tagb and M and Groenewegen (2006) include alu¬ 
minum oxide dust along with the more standard silicate 
in oxygen-rich AGB stars. In future work we will explore 
the effect of different grain species on the mid-IR spectra 
of early-type galaxies. 
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Figure 15. Same as the upper-right panel in Figure [T^ but with 
two models with only Tdiff included. The dashed represents the the 
standard diffuse dust model with Tdiff = 0.001. The dot-dashed line 
represents an otherwise standard model but with the Qpah param¬ 
eter set at the highest allowable value, 10%, and rdiff = 0.003. The 
diffuse dust values were chosen to be as large as possible without 
exceeding the constraints placed 70fim and 160/rm upper limits. 
It shows that diffuse dust does not fit the 6-12/im features seen in 
galaxies like NGC 2300 i.e. early-type galaxies with dominant old 
age stellar populations. 


4.3. Comparison to Low-Metallicity Galaxies 

Low-metallidty galaxies are another class of objects in 
which circumstellar dust may significantly contribute to 
the integrated SEP owing to their g enerally low diffuse 
dust content. Johnson et al. (20131 compared observa¬ 
tions for a subset ot AlNCSl' galaxies with SPS models. 
They found and discussed the possible causes of a 0.2 
dex under-prediction by the models co mpared to Spitzer 
data from the L VL survey presented i n Dale et al. (20091 
(see Figure 5 in Johnson et al. 2013). i'hey considered 
the possibility that circumstellar dust around AGB stars, 
not included in the SPS models, may be the source of the 
discrepancy. 

In Figure [T6| w e plot the offset between observed IRAC 
8fj,m and MIPS 24^m magnitudes and predicted magni¬ 
tudes for models with (blue squares) and without (red 
circles) AGB dust as a function of metallicity. For about 
half the galaxies, the meta llicities are well c onstrained 
by nebular emission lines (Berg et al. 2012). For the 


rest of the sample, met allicities are estimated the mass - 
metallicity relationship (Berg et al. 2012 Lee et al. 2006); 
the resulting metallicities that are broad ly consistent 
with those estimated from HST GMDs (Weisz et al. 
2011). For each galaxy, the corresponding model has 
the same SFH and metallicity. The model photometry 
for each galaxy is based on the estimated SFH for e ach 
galaxy derived from GMDs; see Johnson et al. (2013) for 
details. The default models do not contain any diffuse 
dust. 

In both panels we include linear fits to both sets of 
offsets; in general the offset is reduced with the inclusion 



Figure 16. Comparison of the differences between observed and 
model magnitudes for low-metallicity dwarf galaxies from the 
ANGST survey. Results are shown for models both with and with¬ 
out circumstellar dust around AGB stars. The integrated light 
model predictions utilize the CMD-based SFHs for the ANGST 
galaxies. The lines plotted in each panel are linear fits to the 
points. The average difference between observed and model mag¬ 
nitudes is nearly zero for the Sfim flux but only slightly reduces 
the tension at 24/im. 


of AGB dust. For the 8.0fj,m flux, the average difference 
between observations and models is close to zero. How¬ 
ever, for the 24/im flux, while there is some improvement 
in the offset in all the galaxies, a significant discrepancy 
remains even after including AGB dust. 

We see a very strong correlation between metallicity 
and the size of the offsets at 24^m. The origin of the 
correlation is unclear to us. In Figure we compare 
the residuals between the observations and models that 
only include circumstellar AGB dust and models that 
include AGB dust and diffuse dust at 24/iTO (gpAH = 0.05 
and dust attenuation index a = —1.3). We see that 
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Figure 17. Comparison of the residuals between the predicted 
and observed 24fim magnitude for the low-metallicity ANGST 
galaxies as a function of metallicity. The addition of diffuse dust 
= 0.05, red points) helps resolve the residuals that persist af¬ 
ter including only circumstellar AGB dust (blue points), indicating 
that diffuse dust is important in these galaxies. 

by including diffuse dust we are able to better capture 
the 24/xm flux. However, the correlation with metallicity 
remains. This indicates that joint modeling of UV-IR 
SEDs is necessary for these galaxies to fully understand 
their characteristics. 

4.4. Comparison to Star-Forming Galaxies at z^l 

Finally, we consider the effect of circumstellar dust on 
t he SEDs of actively star -forming galaxies. 


Kelson & Holden (20101 suggested that AGB stars may 
be responsible tor a substantial fraction of the mid-IR 
luminosities in galaxies from z = 0 to z = 2. In Figure [T8| 
we show color-color diagrams o f galaxies at z ^ 1 from 


the FIREWORKS catalog, seeJ^^^se^aL 


details (also, see Figure 2 from Kelson & Ho 


( 2^ fo r 


den 


5010). 


In each panel we include three inodels: a model with the 
standard amount of AGB dust (“lx AGB dust”) where 
the amount of AGB dust included in the model is as 
described in §2.2; a model with a much larger amount 
AGB dust (“lOOx AGB dust”) where we multiplied the 
standard dust optical depth by 100, and a model with no 
AGB dust but a modest amount of diffuse dust (rdiff = 
0.05). Each model assumes Z = Zq, tsf = 2.0 Gyr, and 


fraction of the galaxy colors but fails to span the whole 
extent of the color range. Even when boosting the AGB 
dust optical depth by a factor of 100 (far beyond what is 
plausible) we are unable to reproduce the full color range 
of the data. Furthermore, reducing the metallicity of the 
model to Z = O.4Z0 and thereby increasing the number 
of carbon stars (see Figure!^, and hence the MIR flux, 
only modestly improves the nts of the model to the data. 
However, a modest amount of diffuse dust can readily 
explain the full range of mid-IR colors. It therefore seems 
unlikely that AGB dust alone can explain the full range 

of the mid-IR colors. _ 

This issue has also been discussed in IMelbourne fcl 
Boyer (|20I3 ) where they compared the mid-IR colors 
ot AGB stars in the SMG and EMC with the c olors of 
AGB stars adop t ed for the models presented in |Kelson| 
& Holden| (|2010D. |Melbourne & Boyer| (|2013D concluded 
that th e average AGB color adopted by lKelson fc Holden| 
(2010) were significantly redder than the obse rved aver^ 
age AGB colo rs. This is due to the fact that Kelson & 


Holden (2010) used heavily dust enshrouded AGB stars 


as re presentative of more typ ical AG B stars. Further- 


more, 


Kelson & Holdenl (2010) used the Maraston (20051 


stellar population models, which ap pear to assign too 


much weight to th e TP-AGB phase (|Kriek et al. 
Zibetti et al.||20l5 ). 


2010 


'i'he results ot this section are consistent with our find¬ 
ings from §4 .1 as well as the findings of [Melbourne fc| 
Boyer (2013). Actively star-forming galaxies generally 
contain significant quantities of diffuse dust that tends 
to overwhelm the effect of circumstellar dust in the mid- 
IR. 

5. DISCUSSION 

In this work we computed a set of empirically- 
calibrated radiative transfer models for the dusty cir¬ 
cumstellar shells around AGB stars. These models were 
coupled to a stellar population synthesis code in order 
to explore under what circumstances AGB dust might 
impact the integrated SEDs of composite stellar popula¬ 
tions (i.e., star clusters and galaxies). The models indi¬ 
cate that AGB dust will play an important role in shap¬ 
ing observed SEDs of both quiescent and star-forming 
galaxies when 4/rm A < 8p,m. At longer wavelengths, 
the influence of AGB dust will be substantial only when 
the diffuse dust content is low (Ay < 0.1). Th ese conclu¬ 
sions a re in good agreement with the results of jSilva et al.j 


(1998), who included both circumstellar dust around 
AGB stars and diffuse dust in their SED models. 

Modeling the contribution of dusty circumstellar shells 
to the integrated light of stellar populations carries a 
number of uncertainties that are presently difficult to 
quantify. These include the lifetimes, luminosities, mass- 


the models from 

Kelson & Holden 

(2010 

). In the first 

structure, extent, and lifetime of the circumstellar shell. 

panel we see no ditterence among t 

he ditterent models. 

the types and quantity of dust within the shell, and the 


This is expected since only optical colors are plotted and 
no dust will effect the flux at these wavelengths. This 
plot serves as a check that the models we have adopted 
are a reasonable description of the data and that they 
span the full range of optical-optical colors. 

The second panel is an optical-mid-IR color-color plot 
and we see significant differences among the models. The 
standard AGB model is able to reproduce only a small 


optical properties of dust grains. We attempted to cali¬ 
brate and/or mitigate these uncertainties by tuning the 
model to match several key observables include the SEDs 
of oxygen and carbon-rich stars, the morphology of the 
color-magnitude diagram of the LMC and SMC, the lu¬ 
minosity function of stars in the LMC and SMC, and the 
numbers of AGB stars in low metallicity galaxies from 
the ANGST survey. 
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Figure 18. Color-color plots of galaxies at 5 : ~ 1 from|Wuyts et al.| (|2008|l. The data are compared to stellar population models with 
Tgp = 2.0 Gyr. Three models are plotted in each panel: a model with only diffuse dust (solid blue; rjifi = 0.05), a model with the default 
amount of AGB dust (dashed green), and a model with the AGB dust optical depth increased by a factor of 100 (red dotted). Along each 
line the model age varies from 10 Mvr to 13 Gyr. The mo dels without diffuse dust are unable to match the locus of mid-IR colors for the 
bulk of the galaxies (compare with [Kelson h Holden|2010||. 


The luminosity functions proved to be a strong test 
of the AGB lifetimes in the Padova stellar models. The 
default Padova models produce too many stars even at 
4.5/1771 where dust emission plays a minor role, suggest¬ 
ing that the circumstellar dust models are not the source 
of the discrepancy. After adjusting the lifetimes down¬ 
ward (by lowering their weight in the synthesis) we found 
good agreement for both the d.S/im and S.O^m LFs for 
the LMC. However, the model SMC LF still substan¬ 
tially overpredicts the observed counts by factors of 3 —5. 
Moreover, our recalibration of the Padova models pro¬ 
duces relatively good agreement with AGB counts of low 
metallicity ANGST galaxies, with a typical overpredic - 
tion of only 1.5 (in agreement with Girardi et al.||2010). 
This indicates that the Padova isochrones are doing well 
both at low metallicities probed by the ANGST galax¬ 
ies ( Z < 0.15Zq) and LM C metallicity, but odd ly not 
at SMC metallicity. The Girardi et al. (20101) dusty 
isochrones show a similar overprediction of the observed 
SMC luminosity function. This suggests that there is 
a problem with either the SMC IRAC LF, our adopted 
star formation history for the SMC, or perhaps there is 
some unknown issue with the models unique to SMC-like 
metallicities and ages. 

There is a well-known degeneracy between stellar age 
and metallicity when modeling optical data that is typi¬ 
cally only broken when the strengths of b oth balmer and 
metal lines can be accurately measured. [Bressan et al. 


(1998) argued that circumstellar dust emission around 


AUB stars might also help to break this degeneracy, since 
the contribution of AGB stars to the bolometric luminos¬ 
ity varies strongly with time. We confirm the potential 
power of mid-IR emission as an age diagnostic by fit¬ 
ting ISO spectroscopy of NGC 1404, a galaxy with no 
on-going star formation. The models strongly prefer an 


age of 7.5 Gyr and are capable of discriminating between 
even 1 — 2 Gyr differences in age at old times. Remark¬ 
ably, the age derived in this way agrees very well with the 
analysis of the optical spectrum for this galaxy. The key 
limiting factor to deriving ages based on mid-IR data 
is systematic uncertainties in the modeling. Reducing 
these systematic uncertainties with targeted observations 
of AGB stars in the Local Group would enable firm ab¬ 
solute ages to be obtained in the mid-IR, even for old 
stellar populations. 

With the IRAG and MIPS 24/im p hotometry and 
Syitze r-IRS spectroscopy present ed by [Martini et al.| 
(2013) and Rampazzo et al. (2013), respectively, we con- 
structed a small sample of quiescent, early-type galax¬ 
ies lacking evidence for cold dust. The models presented 
here provide an reasonable fit to the data. This indicates 
that the majority of the mid-IR emission in these galax¬ 
ies is_dominated Iwcircumstellar dust, as also argued 
by Martini et al. (2013). In detail, the models slightly 
over predict the 10/im silicate feature, which may be a 
reflection of our overly simplistic treatment of the grain 
s pecies in the circumst ellar shells. 

Johnson et al. (20I3|) used CMD-based star formation 
histories to make predictions for the integrated light of 
nearby dwarf galaxies. Comparison of these predictions 
to observed total fluxes revealed good agreement in the 
FUV—NIR, but an offset was seen at 8y,m. We repeated 
this test with our updated models and found that the 
offset is reduced from 0.33 dex to 0.13 dex when circum¬ 
stellar dust is included. We also compared the models to 
the observed 24/7777 flux and found that the offset without 
circumstellar dust, which is greater than a factor of 30, is 
only slightly reduced when circumstellar dust is included. 
It is likely that low levels of diffuse dust in these systems 
is providing the bulk of the observed 24/7777 emission. 
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There has been much discussion in the literature re¬ 
garding the sources responsible for the mid-IR emission 
in typical star-forming galaxies. Possibilities include dust 
heated by young stellar populations, dust heated by old 
stellar populations, and circum stellar dust around AGB 


stars. Kelson & Holden (2010) suggested that the latter 
could play a large role in the mid-lR emission of star¬ 
forming galaxies. This possibility is not borne out in our 
models. Instead, circumstellar dust provides a rather 
small contribution to the mid-IR flux for galaxies with 
typical amounts of diffuse dust (Ay > 0.1). This result 
does not preclude the possibility that evolved stars, in¬ 
cluding AGB stars, contribute substantially to th e heat¬ 


ing of the diffuse dust. Indeed, Salim et al. (20091 argued 
for this based on observations of st ar-forming gaL 


1 (see also Utomo et al. 


The launch of JWST wil 


axies at 

20U\. 

turn the spotlight on the 
^ 1 — SO/rm SEDs of stellar populations across time. 
In this wavelength interval a variety of distinct physical 
processes contribute to the light, including PAH emis¬ 
sion and circumstellar dust. Models such as the one pre¬ 
sented in this paper will be necessary to interpret this 
data. In addition, JWST will provide stringent tests of 
these models by observing nearby AGB stars in detail 
and entire populations of such stars in nearby galaxies. 
In the process we hope to learn not only about the evo¬ 
lutionary history of galaxies but also about the physics 
shaping the final stages in the lives of AGB stars. 

6 . SUMMARY 

We now summarize our main results. 


• In agreement with previous work, our models in¬ 
dicate that the dusty circumstellar shells around 
AGB stars can have a factor of ~ 10 effect on 
the integrated flux of star clusters and galaxies at 
A > AfjLm. The effect is especially large at inter¬ 
mediate ages ^ 0.1 — 3 Gyr and when the presence 
of diffuse dust is minimal. The emission associated 
with circumstellar dust shells may be important for 
interpreting the SEDs of quiescent galaxies at high 
redshift, where the average stellar age is relatively 
young. 

• The presence of diffuse dust in a galaxy sharply 
reduces the influence of circumstellar dust emis¬ 
sion at > S/rm. At 24^m, circumstellar dust con¬ 
tributes only ~ 10% of the flux for star-forming 
galaxies with Ay O.I. At shorter wavelengths 
(i.e.,< lO/rm), where under normal conditions dif¬ 
fuse dust does not contribute substantially, circum¬ 
stellar dust may still play an important role. 

• In early-type galaxies lacking signatures of cold 
diffuse dust, the mid-IR data is far in excess of 
a blackbody extrapolation at shorter wavelengths. 
Models that include circumstellar dust shells accu¬ 
rately reproduce the mid-IR photometry and spec¬ 
tra of early-type galaxies. The model MIR flux is 
sensitive to age, even at old ages. This suggests 
that this wavelength range can be a useful age di¬ 
agnostic once circumstellar dust models have been 
accurately calibrated. The lO/rm silicate feature 
is qualitatively well matched but there are some 
quantitative differences at the ~ 10% level. 


We would like to thank Martha Boyer, Martin Groe- 
newegen, Dan Kelson, Ivo Labbe, Sundar Srinivasan, 
Benjamin Sargent, and Roberto Rampazzo for gener¬ 
ously sharing their data, and Paul Martini, Dan Kel¬ 
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on this topic. We also thank the anonymous referee for 
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